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Abstract — Haze is generally phenomenon that an 
atmospheric stuff that greatly degrades the visibility of 
outdoor scenes. This phenomenon occurs due to the 
atmosphere constituent parts that absorb and scatter the light. 
We use the spread map of image in terms of intensity of image 
except a continuous region which has no edge information. 
The scheme performs a per-pixel manipulation, which is 
direct to implement, and then apply the Directed filter to 
improve the image quality. The simulation results 
demonstrate the comparative results and even better than the 
more complex state-of-the-art techniques, having the 
advantage of being appropriate for real-time applications. 
Keywords — Image Dehazing, Blind Dehazing, Direct Image 
filter 

I. Introduction 

Haze is an annoying factor when it shows up in the image since 
it causes poor visibility. This is the major problem of some 
applications in the field of computer vision, such as 
surveillance, object recognition, etc. In order to obtain the clear 
images, haze removal is inevitable. Fog, mist and some other 
particles that degrade the scene image are the results of 
atmospheric absorption and light scattering. The radiance 
achieved to camera along the sightline is decreased due to 
atmospheric light and it is replaced by previously scattered 
light, which is called the air light [1]. This degradation will 
cause the image to lose contrast and color correctness. 
Furthermore, the airlight which affect the image depends on the 
depth of the scene. This information is commonly used for 
dehazing problems. We also adopt this clue to solve the haze 
removal problem. Image haze removal has gotten a growing 
interest recently. More and more methods are introduced in the 
past three years [1, 2]. Nevertheless, dehazing is a challenging 
topic since the haze is dependent on the unknown depth 
information. Often, the images of outdoor scenes are degraded 
by bad weather conditions. In such cases, atmospheric 
phenomena like haze and fog degrade significantly the 
visibility of the captured scene. Since the aerosol is misted by 
additional particles, the reflected light is scattered and as a 
result, distant objects and parts of the scene are less visible, 
which is characterized by reduced contrast and faded colors. 
Restoration of images taken in these specific conditions has 
caught increasing attention in the last years. This task is 
important in several outdoor applications such as remote 
sensing, intelligent vehicles, object recognition and 
surveillance. In remote sensing systems, the recorded bands of 
reflected light are processed [1], [2] in order to restore the 
outputs. Multi-image techniques [3] solve the image dehazing 
problem by processing several input images that have been 
taken in different atmospheric conditions. Another alternative 
[4] is to assume that an approximated 3D geometrical model of 



the scene is given. In this paper of Treibitz and Schechner [5] 
different angles of polarized filters are used to estimate the 
haze effects. A more challenging problem is when only a 
single degraded image is available. Solutions for such cases 
have been introduced only recently [6]-[9]. In this paper we 
introduce an alternative single-image based strategy that is 
able to accurately dehaze images using only the original 
degraded information. An extended abstract of the core idea 
has been recently introduced by the authors in [10]. Our 
technique has some similarities with the previous approaches 
of Tan [7] and Tarel and Hautiere [9], which enhance the 
visibility in such outdoor images by manipulating their 
contrast. However, in contrast to existing techniques, we built 
our approach on universal dehazing with directed filter. We 
are the first to demonstrate the utility and effectiveness of a 
fusion-based technique for dehazing on a single degraded 
image then we built the universal image dehazing model with 
directed filter. In this work, our goal is to develop a simple 
therefore; all the universal dehaze processing steps are 
designed in order to support these important features. The 
main concept behind universal dehaze based technique is that 
two input images from the original input with the aim of 
recovering the visibility for each region of the scene in at 
least one of them [7]. Additionally, the universal dehaze 
image enhancement technique estimates for each pixel the 
desirable perceptual based qualities (called weight maps) that 
control the contribution of each input to the final result. In 
order to derive the images that fulfil the visibility assumptions 
(good visibility for each region in at least one of the inputs) 
required for the fusion process, we analyse the optimal model 
for this type of degradation. 

II. Haze Detection by Universal Dehazing Method 

The human eyes sense more brightness than color. Therefore 
the atmospheric light valuation uses and create a transmission 
map in the YC- 0 C r color channels. The atmospheric light is 
predictable from the most opaque pixel. 
In the dark channel prior the existing procedure picks up the 
top 0.1% brightest pixels. Sine an image does not have 
information on the edge of the sky or a wall in the area, the 
mis-estimated value of the atmospheric light results in failure 
of the defogging (dehazing) algorithm. Therefore to represent 
the neighboring pixel's relative depth information we use the 
edge information. We can construct the consistent 
atmospheric light with this relative depth information to 
restrain the edge halation. With this relative depth 
information map we generate the transmission by estimating 
the atmospheric light except a continuous region which has 
no edge information. And the transmission map is given as, 
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Where £ ^restricts the transmission t (x) to a lower bound t D , 
which means that a small amount of fog are preserved in very 
dense fog regions. Color distortion problem may occur in the 
compensation process. To solve this problem, the image 
restored by color correction using statistical RGB channel 
feature extraction of image. We calculate the RGB channel 
ratio between foggy and defogged images for color correction 
with weighted image. The RGB channel ratio is defined as, 
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Where R represents the defogged image an O the foggy image, 
As a result, we can obtain the color-corrected image using color 
matching of RGB channels of restored image, such as 



1 = 
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Where / represents the color-corrected image, and k the number 
of pixel 

A. Directed filter Image modelling for Haze 
Extraction 

The observed brightness of a capture image in the presence of 
haze can be modelled based on the atmospheric optics [6, 7,11] 
via 

/(*) =J{ X }t(x) + A(l -£(*)} (3) 
Where, I(x) is the observed haze image, J(x) is scene irradiance 
(the clear haze-free image), A is the airlight that represents the 
ambient light in the atmosphere. t(x)c[0, 1] is the transmission 
of the light reflected by the object, which indicates the depth 
information of the scene objects directly. J(x)t(x)on the right 
hand side is called direct attenuation, which describes the scene 
radiance and its decay in the medium. The second term A(l- 
t(x))is the atmospheric veil (atmospheric scattering light), 
which causes fuzzy, color shift, and distortion in the scene. The 
goal of haze removal is to recover J(x) , A and t(x"j from 

Image Dehazing 

In this section, we will describe in detail. The rough down- 
sampled transmission and the air-light are estimated firstly, 
then the transmission is smoothed and up sampled using a 
directed filter, and finally the haze-free image is restored. 
A. Extract the Transmission 

The core of haze removal for an image is to estimate the 
airlight and transmission map. Assuming the airlight is already 
known, to recover the haze free image, the transmission map 
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should be extracted first. He et al. [8] found that the minimum 
intensity in the non-sky patches on haze free outdoor images 
should have a very low value, which is called dark channel 
prior. Formally, for an image /, the dark channel value of a 
pixel x is defined as: 

/^™(x) = miB^^jCmiBy^^ (j c Cy))) 
where,/ c is a color channel of / ; Q(x) is a patch around x. By 
assuming the transmission in a local patch is constant and 
taking the min operation to both the patch and three color 
channels, the haze imaging model in (4) can be transformed 
as: 

mil W^b) ( miT V^a W (^r))+ O^W) ( 4 ) 
where, £(x) is the patch transmission. Since A is always 
positive and the dark channel value of a haze-free image / 
tends to be zero according to the dark channel prior, we have 

Then the transmission can be exacted simply by: 

F(x) = min^^ (min y£flM (^r)) ( 5 ) 
Although the dark channel prior is not a good prior for the sky 
regions, fortunately, both sky regions and non-sky regions can 
be well handled by (5) since the sky is infinitely distant and 
its transmission is indeed close to zero. In practice, the 
atmosphere is not absolutely free of any particle even in clear 
weather. Therefore, a constant parameter co(0<co<l) is 
introduced into(8) to keep a small amount of haze for the 
distant objects: 

f (x) = 1- co miiw^ (min^u (^r)) (6) 
The estimated transmission maps using (6) is reasonable. The 
main problems are some halos and block artifacts. This is 
because the transmission is not always constant in a patch. 
Several techniques were proposed to refine the transmission 
map, such as soft matting and directed joint bilateral filter. 
These techniques were applied on the transmission maps of 
the original foggy images and usually several operations 
should be used to achieve a good result, which could be 
computational intensive. For image haze removal, the time 
complexity is a critical problem that needs to be addressed. 
High time complexity of dehazing may make the algorithm 
impracticable. 
B. Refine the Transmission 

To improve the efficiency, in the present implementation, the 
transmission map is obtained from a down-sampled minimum 
channel image. Then, it is refined and up-sampled by using 
directed filter, which can be explicitly expressed by [11]: 

t^ljW^ij (7) 

Wij (/*) = j^S^y^Cl + ( 8 ) 
Where, J^is the guidance image; u-^and a-* are the mean and 
variance of /^inw^ ; Iwl is the number of pixels in w^.f is a 
regularization parameter. The refined operation on a down- 
sampled minimum channel image leads to a low time 
complexity and helps to reduce halos and block artifacts. Joint 
up sampling using directed filter is applied to obtain the full 
transmission map. The directed filter is reported to be a fast 
and non-approximate linear-time algorithm, which can 
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perform as an edge preserving, smoothing operator like the 
bilateral filter, but does not suffer from the gradient reversal 
artifacts. Moreover, the directed filter has an O(N') time (in the 
number of pixels N) exact algorithm for both gray-scale and 
color images. 

B. Work Flow of the Algorithm 

Image acquisition: Take the input image as we perform all 
possible operation on it. Figure 1 showing all the procedure 
related to haze removing algorithm. 
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possible configurations. An alternative method for 
approximation to ML estimation is maximum pseudo 
likelihood (MPL) estimation, which was proposed by Besag 
[15]. The MPL estimation method is a suboptimal method, 
which is given by 

& = arg miti c T s V c (H(s~) , a) - InZ^ . (11) 
The differences among the fused results are usually difficult 
to be measured only based on observation, particularly when 
the fused images are multiband. Objective and quantitative 
analysis can benefit to a comprehensive evaluation. Various 
image quality indices have been developed for the purpose of 
image fusion [12]— [13]. Some of these indices validate the 
spatial resolution, while others focus on the spectral 
properties of the obtained fused result. In this paper, we 
employ three such indices. 

1) SNR: The SNR in decibels, as shown in (12), is a direct 
index to compare the fused image to the reference one 
[16] .For multiband images, it can be calculated band-by-band 
and also globally averaged SNR 

SNR (Z, Z) = lOlogio s ff )g ( 12 > 
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Figure 1 Flow chart and Algorithm for Haze Removal 

IB. Result Analysis and table 

The algorithm proposed here will remove haze from an image 
surface without prior knowledge of the haze location upon that 
surface. The proposed method is based on determining the 
illumination profile of the image surface. There are some 
Performance Parameters used for algorithm comparisons values 
of these evaluation metrics should be high. 
Modelling the Markov pdf parametrically involves the data 
driven optimal estimation of the parameters associated with the 
potential functions Vc. The model parameters must be 
estimated for each data set as part of the image processing 
algorithm. In our algorithms, the noise variance ol in (6) and 
the parameter a in the coefficient MRF pdf in (7) are unknown. 
Thus, we need to estimate these parameters in our algorithms. 
Because we assume that the noise in the fusion model is a 
Gaussian noise, it is straightforward to estimate the noise 
variance by the maximum likelihood (ML) criterion. It is given 
by 

ff2 = ^LGVtfL*f tt-fW (9) 

The direct ML estimation of the parameters associated with the 
pdf of H is known to be a difficult problem [32]. The ML 
estimate of a is 

a = irg mEti c V C (H, a) — lnZ H (10) 
The potential function Vc(H, a] can be simply computed. 
However, the normalization term ZH involves a summation 
over all possible configurations of H, which is practically 
impossible due to the large computation time. Note that, for 
two source images with size 300 *300, H has a total of 490000 



2) Universal Image Quality Index (UIQI): A UIQI [14] has 
been widely used for image similarity evaluation and was also 
applied to validate fusion techniques [13]. UIQI of two 
images (A and B) is defined as 



Q = 



(13) 



This quality index models any distortion as a combination of 
three different factors: loss of correlation, luminance 
distortion, and contrast distortion. The dynamic range of Q is 
[-1, 1], and the best value 1 is obtained if A = B. When 
applying this index to a multiband image, it is applied band- 
by-band and averaged over all bands. [16]. 
3). Performance of the image compression coding, it is 
necessary to define a measurement that can estimate the 
difference between the original image and the decoded image. 
Most image compression systems are designed to minimize 
the MSE and maximize the PSNR. 



MSE = 



PSNR =201og 1D g| 



(14) 



(15) 



This profile is then used to remove the haze. It is 
implemented using MATLAB 7.9.0 (R2009b) on i-5 
processor with 4-GB RAM. The simulations have been tested 
on aerial images in figure 2 and 3. Figure 2 shows the 
Original Image of Building image of haze Removed Image. 
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Figure: 2 (a) Original Image of building, (b) Dehazed Image by 
using Multi scale Fusion (c) Dehazed Image by using Universal 
Dehazing (d) Dehazed Image After Directed filter. 
Table 1 Comparison parameters for image of building 



Method 


Variance 


Mean 


SNR 


UIQI 


Multiscale 
Fusion 


0.1252 


0.5415 


6.4288 


2.8219 


Universal 
Dehazing 


0.1446 


0.4542 


6.5875 


2.9817 
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Figure: 3 (a) Original Image of field, (b) Dehazed Image by 

using Multi scale Fusion (c) Dehazed Image by using 
Universal Dehazing (d) Dehazed Image After Directed filter 

Table 2 Comparison parameters for field image 



METHOD 


VARIANCE 


MEAN 


SNR 


UIQI 


Multiscale 
fusion 


0.1039 


0.5315 


5.7040 


2.7719 


Universal 
dehazing 


0.1241 


0.4352 


7.5421 


2.5927 



IV. Conclusion and Future Scope 

Here in this paper, a fast and effective scheme for image and 
video dehazing is proposed. Using a newly presented haze 
removal for a single image without using any additional 
information is framed as a specific filtering problem and an 
improved filtering scheme is planned based on directed filter. 
In this offered algorithm, the airlight and the down-sampled 
transmission can be estimated and extracted easily. Then 
using a directed filter, the transmission can be further refined 
and up-sampled. Results demonstrate the presented method 
abilities to remove the haze layer and achieve real-time 
performances. It is believed that many applications, such as 
outdoor surveillance systems, intelligent vehicle systems, 
remote sensing systems, graphics editors, etc, could benefit 
from the proposed method. 
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